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Critical Dynamics of Three-Dimensional Spin Systems with Long-Range Interactions 
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A field-theoretic description of critical behavior of Ising systems with long-range interactions is 
obtained by using the PadeBorel summation technique in the two-loop approximation directly in 
the three-dimensional space. It is shown that long-range interactions affect the relaxation time of 
the system. 
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It was shown in [1] that effects due to long-range in- 
teraction are essential for the critical behavior of Ising 
systems. The renormalization-group approach to spin 
systems with long-range interactions developed in [2] di- 
rectly in the three-dimensional space made it possible to 
calculate the static critical exponents in the two-loop ap- 
proximation. However, analogous calculations of critical 
dynamics have never been performed for these systems. 

In this paper, a field-theoretic description of critical 
behavior of homogeneous systems with long-range inter- 
actions is developed directly for D — 3 in the two-loop 
approximation. The model under analysis is the classi- 
cal spin system with exchange integral depending on the 
distance between spins. The corresponding Hamiltonian 
is 



H 



)SiSj, 



(1) 



where Si is a spin variable and J{\rirj\) is the exchange 
integral. This model is thcrmodynamically equivalent to 
the 0(n)-symmetric GinzburgLandau Wilson model de- 
fined by the effective Hamiltonian 



H 



d D q 



{\(T + q a )<p 2 + u Q p 4 }, 



(2) 



where <p is a fluctuating order parameter, D is the space 
dimension, r ~ \T — T c \ ( T c is the critical tempera- 
ture), and wo is a positive constant. Critical behavior 
strongly depends on the parameter a characterizing the 
interaction as a function of distance. It was shown in [3] 
that long-range interaction is essential when < a < 2, 
whereas systems with a — 2 exhibit critical behavior 
characteristic of short-range interactions. For this rea- 
son, the analysis that follows is restricted to the case of 
< a < 2. 

Relaxational dynamics of spin systems near the crit- 
ical temperature can be described by a Langevin-type 
equation for the order parameter: 



dip _ _ x SH 
dt Sip 



h i] + A h, 



(3) 
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where Ao is a kinetic coefficient, rj(x, t) is a gaussian ran- 
dom force (representing the effect of a heat reservoir) 
defined by the probability distribution 



P n = A v exp 



-(4A0)" 1 / d d xdtr) 2 {x,t) 



(4) 



with a normalization factor A n , and h(t) is an external 
field thcrmodynamically conjugate to the order param- 
eter. The temporal correlation function G[x, t) of the 
order-parameter field can be found by solving Eq. (3) 
with H[ip] given by (2) for <p[?7,h] averaging the result 
over P v , and retaining the component linear in h(0) 



G{x,t) 



<Sh(0) 



[vV>0M))]U=o, 



where 



{( V (x,t))} = B- 1 / D{r,}<p(x,t)P v 



B = J D{ V }P V . 



(5) 



(6) 



(7) 



When applying the standard renormalization-group 
procedure to this dynamical model, one has to deal with 
substantial difficulties. However, it was shown in [4] that 
the model of critical dynamics in homogeneous systems 
without long-range interaction based on a Langevin-type 
equation is equivalent to that described by the standard 
Lagrangian [5] 



L = d d x dt 



up 



-idip 
dt 



Sip 



(8) 



where ip* denotes an auxiliary field. The corresponding 
correlation function G(x, t) of the order parameter is de- 
fined for a homogeneous system as 

G(x,t) = (ip(0,0)ip(x,t)) (9) 

IT 1 / D{<p}D{<p*}<p{0,Q)<p(x,t)exp(-L[<p,<p*]), 



where 



n = J D{<p}D{<p*}ejq>(-L[<p, ip*}). 



(10) 
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Instead of dealing with the correlation function, it is rea- 
sonable to invoke the Feynman diagram technique and 
represent the corresponding vertex in the two-loop ap- 
proximation as 



r (2) (fc, w; T , u , An) = T + k a - — - 96u 2 D , 



d D qd D p 



(i + i?r)(i + w) 

1 



(11) 



(3 + W + W + W + o\ c 



VA) 



The next step in the field-theoretic approach is the cal- 
culation of the scaling functions (3, 7 T , 7^ and 7a in the 
renormalization-group differential equation for vertexes 



— +13— 
djj, du 



lr-r- 



d 



7aA 



d 



m 



xT^(k,w,T,u,X,n) = 0, 



-1<p 



(12) 



where the scaling parameter /i is introduced to change to 
dimensionless variables. 

Further analysis requires the use of the function (3 and 
the dynamic scaling function 7a. 

An expression for (3 in the two-loop approximation 
were obtained in [2]: 



[3 = -(4 - D) 1 - 36uJ + 1728 (2 Ji - J 2 - -G)u 2 ) 

d D qd D p 



Jl 



(1 + |<fr ) 2 (1 + \p\°)(l + \q 2 +p 2 + 2pq}°/ 2 ) ' 
f d D q 

J (i + WY 

G = 



d 



d D qd D p 



d\k\° J [l + \q 2 + k 2 + 2k^/ 2 ){\ + \^) 
1 



[l + \q 2 +p 2 + 2pq}°/ 2 ) 



The function 7a calculated in the two-loop approxima- 
tion is 



7a = (2a - D)2(D' - G)u 2 , 
dDo 



(13) 



D' 



d(-iiu/X) 



Defining the effective interaction vertex v — u ■ Jq one 
obtains the following expressions for (3 and 7a: 



[3 = -(2a -D) 
7A = (2cr - D)96(D - G)v 2 



1 - 36w + 1728 (2 Ji - 1 - -g)« 2 ) 



(14) 



Ji 



G 



D' 



Jl — —f2 G — —[2 D — — 2. 

J Q J Q Jq 



This redefinition is meaningful for a < D/2. In this case, 
Jo, Ji, G and D' are divergent functions. Introducing the 



cutoff parameter A, we obtain finite expressions for the 
ratios 



J 2 , 



II 

JO JO 
(l + \q 2 +p 2 + 2pqr)) 



d u qd u p/{{\ + w r(i + W) 

A 



G_ 

I 2 
J o 



d 



A r A 



d u q /{i + \qry 



d u qd u p/{{\ + \PT) 



(d\k\°)Jo Jo 
(1 + \q 2 + k 2 + 2kqX )(1 + \q 2 + p 2 + 2pq\ a )) 
2 



/ d D q/{l + \q\°f 
Jo 



(15) 



D' _ 3 
J 2 "4 



rfVV((i + kT)(i + bT) 

A 



(3 + i<zr + bT + ip+9r) 



as A — > 00. 

The integrals are performed numerically. For a < D/2, 
a sequence of Ji/Jq and G/Jq corresponding to various 
values of A is calculated and extrapolated to infinity. 

Critical behavior is completely determined by the sta- 
ble fixed points of the renormalization group (RG) trans- 
formation. These points can be found from the condition 



(3{v*) = 0. 



(16) 



The effective interaction vertexes were evaluated at the 
stable fixed points of the RG transformation in [2] . 

The dynamic critical exponent z characterizing critical 
slowing-down of relaxation is determined by substituting 
the effective charges at a fixed point into the scaling func- 
tion : 



2 = 2 + 7A . 



(17) 



The table shows the stable fixed points of the RG trans- 
formation and the values of the dynamic critical expo- 
nent for 1.5 < a < 1.9. When < a < 1.5, the only fixed 
point is the unstable gaussian one, v* = 0. 

A comparison of the present results with the value 
z = 2.017 of the dynamic critical exponent for three- 
dimensional systems with short-range interactions ob- 
tained in [6] demonstrates the essential role played by 
long-range interactions in critical dynamics of spin sys- 
tems. In particular, the system's relaxation time in- 
creases according to the scaling t ~ \T — T c \ vz , where 
v is the critical exponent characterizing the increase in 
the correlation radius near a critical point. In three- 
dimensional systems with long-range interaction, both 
critical dynamics and static behavior become increas- 
ingly gaussian as the long-range interaction parameter 
a decreases. When a < 1.8, critical behavior is virtually 
gaussian. 

The work is supported by Russian Foundation for Ba- 
sic Research N 04-02-16002. 
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Tabic 1. 



a v* z 

1.5 0.015151 2.000072 

1.6 0.015974 2.000180 

1.7 0.020485 2.000777 

1.8 0.023230 2.001529 

1.9 0.042067 2.006628 
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